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Abstract 

We propose an experiment in which the phonon excitation of ion(s) in a trap, with a trap 
frequency exponentiahy modulated at rate k, exhibits a thermal spectrum with an "Unruh" tem- 
perature given by T = Hk. We discuss the similarities of this experiment to the usual Unruh 
effect for quantum fields and uniformly accelerated detectors. We demonstrate a new Unruh effect 
for detectors that respond to anti-normally ordered moments using the ion's first blue sideband 
transition. 
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I. INTRODUCTION 



It has been known for many decades that an accelerated detector, moving in a quantum 
field prepared in the ground state of the field modes for an inertial frame, will become 
excited|l|,|2,|3[. In the case of constant acceleration, a, the frequency response of the detector 
is completely equivalent to the response of an inertial detector in a thermally excited field 
with a temperature T given by T = ha/2T{c. We call this Unruh-Davies radiation. Quite 
clearly, such an effect would be very difficult to see given current technologies. In this paper 
we suggest an analogous system, based on detecting phonons of the vibrational modes of 



cold trapped ionsj^, |5[. In many ways this parallels a theme, pioneered by Unruh, of sonic 
equivalents for quantum fields in curved space time ^ . 

Our analogy is based on an alternative view of Unruh-Davies radiation in terms of the 
time dependent red shift seen by an accelerated observer j.?]. By controlling the trapping 
potentials of trapped ions it is possible to modulate the normal mode frequencies so that 
they have the same time dependent phase as red-shifted frequencies seen by a constantly 
accelerated observer. Suppose now that the ions are prepared (using laser cooling) in the 
ground state of the normal modes of vibration of the time independent trap. If a suitable 
detector of the vibrational quanta for the trapped ions was available, they could be used 
to detect excitations out of the ground states of vibrational motion due to the frequency 
modulation. This would be analogous to the response of an accelerated detector to a scalar 
field prepared in a Minkowski ground state. Fortunately the ions themselves can be made 
to respond as phonon detectors. 

In ion trap implementations of quantum information processing, a laser is used to couple 
the vibrational motion of a trapped ion to an electronic transition between states which we 

n 

denote ^ e [8|]. Furthermore it is possible to readout one or the other of these internal 
electronic states, say using a laser (the readout laser) to drive a cycling transition between 
g and another electronic state, thereby producing fiuorescence conditional on whether the 
ion was in state g. Such measurements are highly efficient and closely approximated by 
a perfect projective measurement of the internal electronic state. In effect this scenario 
defines a phonon detector that may be turned on and off at will. To be more specific, we can 
implement various kinds of phonon detectors by carefully tuning the laser frequency ujl to 
one of the vibrational sidebands of the ion. This enables one to realize rather unconventional 
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phonon detectors that respond to antinormally ordered moments (blue sideband) as well as 
the more conventional normally ordered moments (red sideband) as we explain in more 
detail below. 

II. TRAPPED ION MODEL. 

The interaction Hamiltonian describing the coupling of the internal and vibrational de- 
grees of freedom of the m'th ion in a linear array of ions in a trap (in the interaction picture, 
and Lamb-Dicke limit) can be written as 

= hnok cos eqm{t)cT^{t) (1) 

where 

a..(t) = e-^^V_+h.c., (2) 

and Qq is the (scaled) Rabi frequency, A = ua—ujl is the detuning between atomic resonance 
and the laser, k = ul/c is the wavevector, 6 is the angle the laser beam makes with the 
longitudinal axis for the linear ion chain, qm{t) is the quantized local displacement of the mth 
ion about its equilibrium position and cr_ = \g){e\ is the atomic lowering operator between 
the upper and lower atomic states |e) and \g) separated by frequency ua- We assume the 
equilibrium position of the mth ion is located at a node of the standing wave laser beam, 
we have used the rotating wave approximation to describe the interaction between the laser 
and the ion. In the Lamb-Dicke limit, terms of order q'^{t) have been neglected since the 
displacement of the ion is much less than the wavelength of light. In an experiment the 
e ^ g is a quadrupole transition and is driven by a Raman process. Thus Qq is an effective 
Rabi frequency for this process. 

Equation (jT)) could also be regarded as a discretised representation for the interaction of 
a scalar field, (pix), and a local detector, with transition frequency A, where (f){x = ml) = 
for some discretisation length /. In such an interpretation the interaction Eq.(^is equivalent 
to the Unruh model of a particle detector 0, 0], with only two internal energy levels. 

Note that the detector can be turned on and off through the dependance on the external 
laser field in fig, a somewhat unusual feature for field quanta detectors. Another unusual 
feature of this detector is that the transition frequency of the detector, A = tu^ — u;^ can 
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be varied by tuning the external laser. Conventional detectors would have a fixed transition 
frequency. This latter feature will enable us to define different kinds of phonon detectors. 

The local displacement of the mth ion, qm{t), can be expanded in terms of creation and 
annihilation operators for global normal modes (phonons) of the A^-ion system Qp by ^ 

N 

p=i 

— n — ^ 

^T.^'lH-pe-^' (3) 
p=i 

where the coupling constant is defined by 

^ (4) 

In the above Up are the normal mode trap frequencies given by Up = ^JJipU in terms of 
the bare trap frequency v and the eigenvalues [ip. For the center of mass mode s^m = 1, 
ui = u and for the breathing mode z/2 = y/3iy, Sm = VN {^^^i Um / v^, where in 

the later Um are the components of the breathing mode eigenvector bm = Uml \jYlim=i ""m 
normalized to unity. Further normal mode parameters for up to = 10 ions are computed 
numerically in James 



Defining the Lamb-Dicke parameter as ?7 = a/ hk"^ cos^ 9/2mv we can write our Hamilto- 
nian in its final form 

= xUtM) (5) 

where we have defined the operator 



N 



* p=i 

and X = ^o'?- This Hamiltonian describes the coupling between the two level electronic 
transition (the detector) and the vibrational degrees of freedom whenever the external laser 
is turned on {Qq ^ 0). Note that we do not make any assumptions at this stage about the 
relative size of A and Vp. We wish to keep A as a free parameter which may be varied to define 
different kinds of phonon detectors. In analogy with the Unruh detector model, the field 
operator qrrSf) represents the scalar field at the position x = ml and time t, qm(t) = (j){x,t). 
In physical terms this is the displacement of the mth ion as a function of time. 
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It is worth noting an important difference between this model and the usual treatment of 
a particle detector. In the case of a usual detector the frequency term, A, would be strictly 
positive thus defining the positive and negative frequency components of the dipole a±. In 
Eq.Q, the parameter A can be positive or negative so we cannot simply refer to positive 
or negative frequency components in absolute terms. However the operators a± will retain 
their usual definition as raising and lowering operators. 

In the case that A > 0, the laser is detuned below the atomic transition, which we refer to 
as red detuning. We can resonantly excite so called red sidband transitions when A ^ nvp. 
Near such a resonance (n = 1) we can make the rotating wave approximation and describe 
the interaction by the Hamiltonian 



This describes the usual Jaynes-Cummings model of a two level system interacting with a 
bosonic degree of freedom. In physical terms it describes a Raman process in which one 
laser photon and one trap phonon are absorbed to excite the atom (see figure CJ). A phonon 
detector defined this way would respond to the normally ordered moments of the phonon 
field amplitude. 

In the case that A < 0, the laser is detuned above the atomic transition, which we will 
refer to as blue detuning. The resonant term for the fisrt blue sideband is then given by 



Again this is a Raman process in which the atom is excited by the absorption of one laser 
photon and the emission of one phonon (see figure P). Considered as a phonon counter this 
would correspond to a detector that responded to anti-normally ordered moments of the 
phonon field amplitude. 

Using laser cooling techniques it is possible to cool the system very nearly to the ground 
state of the vibrational degrees of freedom. In reality this becomes more difficult as the 
number of ions, and thus normal modes increases. However for our purposes even one ion 
would suffice. In current experiments the cooling is sufficiently efficient to reach the ground 




p=i 





p=i 





^(0)) = |0) = \g) ® |0)vib 



(9) 
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where the initial vibrational state is a tensor product of the ground states of each of the 
normal modes, 

p 

io)vib = n • (10) 

n=l 

We will exponentially chirp the trap frequency up or down such that 

p jj{t) = e^"' (11) 

with K the chirp rate and focus our laser on the first ion (m = 1). For a constant trap 
frequency the phonon annihilation operator satisfies the usual uncoupled mode equation 
ap{t) = —iVpap{t) with solution ap{t) = exp{—ti>pt) ap{0) , Vp = y/Jj^i' which was used in 
Eq.®. For an chirped trap frequency given by Eq. pi|) between an initial time to and final 
time t, the phonon mode now satisfies the following equation of motion, with solution 

^dp{t) = -iupe^^ap{t) ap{t) = exp (z^e'^*") exp (-z^e'^*) ap{to). (12) 

where we consider the chirp-up case first (with the chirp down case discussed subsequently). 

We now suppose that the coupling to the detector is turned on at the same time as 
the frequency modulation is turned on, and turned off at the same time as the frequency 
modulation is turned off. This is a rather different scenario to the usual discussion of the 
Unruh-Davies effect in which the detector is always coupled to the field and continuously 
accelerated, ie the red shift frequency modulation is always on. 

We are interested in the probability Pm{T, to) for the excitation of the m'th ion from the 
ground state to all excited states of detector and field for a detector turned on at t = to and 
turned off at t = T. The detector in our case has only one excited state. We let the excited 
states of the vibrational degree of freedom be represented by a complete orthonormal basis 



i). The total excitation probability is then 10 1 

rT pT 



PUT,to)=x'[ dt' ! t/tV^(*'-*")G(t',t") (13) 

J to J to 

where the field correlation function is defined as 

G{t',t")= vib(0|g„(t')gm(t")|0)vib (14) 

with the field now given by 

N 

^™(^) = ^ E (a.e-^-' - ate^^-') , (15) 



N 1 
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and we have neglected the phase factor arising from the initial time as it does not contribute 
to the correlation function. Substituting this result into Eq.lfT^. 

^ |/,(P)|2 

P».(T,to) = x'E^IW^o)r (16) 



where the integral is 



p=l V 



T 



J,(r,to)= / dte^^'e^'^^^'). (17) 

J to 

With no loss of generality we may now set to = 0. We first consider the case A > 0, which 
is the red-sideband case. Our objective is to calculate the excitation probability for the two 
level system near the red sideband transition which we label P^(T, to)- 

To evaluate this integral we first change to dimensionless time, r = nt, so that 

1 r'^'^ 

Ip{T,to) = - tire*'^V^\ (18) 



K 







where a = A/k, and b = Up/n. Next we define the constant a 

b = e"" (19) 

It then seems sensible to make the new change of variable r' = r — a so that the integral 
becomes 



kT-o 

dre'^^e'^^ (20) 

K 

In an experiment one needs to vary A near the red or blue sideband, so we expect that 
A and Vp are of the same order. We now consider the limit in which Vp/ n << 1 for all the 
normal modes. In this limit a >> 1, the integral does not depend much on a, so we extend 
the lower limit of the integral to infinity. Furthermore we suppose the time over which the 
detector and modulation are on is such that T » Up, in which case the integral does not 
vary much with T. In that case we can extend the upper limit of the integral to infinity, so 
that 

„iaa poo 

I^r^ / rfre^^'e^'^" (21) 

1^ J-oo 

Changing the variable of integration to y = e'^ we obtain an expression for Ip with arbitrary 
limits as 

Ip = / dyy'^-'e'y (22) 

^ Jo 

iaa 

= — r(2a)e-™/2 (23) 

K 
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If we now use the identity |r(ia;)p 



7r/[x sinh(7rx)] for real x ll| we see that 




1 



(24) 



kA e' 



Inserting this into Eq. lfT^ we can write in the suggestive form 



m 



N 



(25) 



p=i 



where we have defined the Unruh temperature as 




(26) 



Eg. ()25j) has the analogous form to a thermal spectrum at temperature T as seen by a uni- 
formly accelerated observer moving through a particle-free inertial vacuum. The "thermal" 
form of the probability is independent of the phonon frequencies Vp since we have taken the 
upper limit t to infinity. The last summation expression is just a numerical factor, which 
can be computed j^, or dropped in the case of a single ion in the trap. 

The analogy with the Unruh effect for a umformly accelerated observer in Minkowski 
space can be seen as follows. As discussed in [7j the chirping of the trap frequencies can 
be considered as arising from the the modification of the usual Minkowski plane wave 
exp(iy9-|-) = exp {i{kz ± u;t)) due to the motion of the accelerated observer. For an observer 
moving at constant velocity Vq in Minkowski space, a Lorentz transformation (LT) of the 
phase ^p± viz ct' = ct cosh r + z' sinh r and z' = z cosh r + ct' sinh r with the rapidity defined 
by tanhr = fo/c, simply transforms it to (p± Lp'^ = {uj/c)e^''' [z! ± ct') = {oj' /c) {z' ± ct') 
which merely produces the constant doppler shifted frequencies in the new frame u' = uoe^^ . 
For an accelerated observer, one has to perform a time dependent LT at each instant 
to the comoving frame that is instantaneously at rest with respect to the accelerated 
observer. The orbit of the accelerated (Rindler) observer as described by an inertial 
Minkowski observer is given by the Rindler transformation ct(r) = (c^/a) sinh(ar/c) and 
z{t) = (c^/a) cosh(ar/c) where a is the constant uniform acceleration. Under this trans- 
formation, the Rindler observer experiences phase of the Minkowski plane wave exp{iip±) 
as transformed to ip± V^±(''") = kz^r) ± cot^r) = (cuc/a) e^""^/^. These are the chirped 
frequencies that appear in Eq. ()17p in exponentially expanding or contracting the trap. The 
Fourier transform of this modified plane wave with respect to the Rindler proper time t can 
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be considered as a measure of the noise spectrum seen by the accelerated observer viz 



SiA) 



j-g±j7r/2^±iAc/a / 



2nc 



Aa e^""!"- - 1 '' 



(27) 

analogous to |/p(c)0, — oo)p as in Eq.(j21I). Here the perceived thermal temperature is defined 
from the only other frequency that can be formed from the Rindler observers accelerated 
motion, a/c. Thus we get the usual Unruh temperature defined by T = Uajl-nc. In our ion 
trap analogy, the role of the acceleration frequency a/c is played by the trap expansion rate 

K. 



0) 



g 



V 



CO, 



g 



A>0 



A<0 



FIG. 1: Energy Level Diagrams. Tuning to (a) red side band, A = loa — (^l = i^p > 0, (b) blue 
side band A = — < 0. The Unruh effect Ea. (|25jl takes place on the red sideband. We obtain an 
anti-normally ordered Unruh effect if we tune to the blue side band. 



It is important to note that in order to obtain the Planck factor (e^'^'^/" — 1)^^ in Eq.( 
indicative of a Bose-Einstein (BE) thermal distribution and the signature of the Unruh effect, 
we made crucial use of a positive detuning A > 0, corresponding to a red sideband detuning 
in Fig.© and Fig.©. This resulted in the factor (i)'^^^ = (e*'^/^)'^/'' = g-^^/^^ that 
appears in Ip in Eq. ()23|) . Dividing the square of this factor into the sinh function appearing 
in the denominator of |/pP, resulting from the term |r(iA//t)p, produces the signature BE 
thermal distribution function. 

If on the other hand, we had instead chosen A = — |A| < corresponding to a negative 



A = 




A>0 





A<0 



FIG. 2: Vibrational Resonances. Red (a) and blue (b) vibrational sidebands corresponding to the 
energy level diagrams in Fig.([T|). 



g7r|A|/2ft I3ividing the square of this term into the sinh function appearing in the denominator 
then produces alternatively the probability for excitation on the blue sideband 



with the same definition of the Unruh temperature as in Eq. (j26p . 

We might label such a distribution an anti-normally ordered Unruh effect since the vi- 
brational excitation from the ground to the excited state takes place by an absorption of a 
photon and an emission of a phonon to the electronic state |e) as depicted in Fig.(P)b, i.e. 
by a term such as apcr_, see Eq.®. As discussed earlier, such a detector responds to the 
anti-normally ordered moments of the phonon field amplitude. Note that as k (T ^ 0) 
in Eq.(|28|) we get a finite contribution to the probability for excitation P^- In the case of 
red sideband detuning, the usual Unruh effect analogy in Eq. lj^ . — > as k — > 0. This 
limit corresponds to a fixed trap frequency z/ for which we get no excitation as described 
above. 

The above k —>■ limiting cases can also be understood as follows. For a constant trap 



detuning to the blue side band Fig.Qb, the previous factor would become (e*'^/^) 




(28) 
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frequency t'p(t) = Vp the relevant integral to compute is 

/oo 
dte'^'e'''^' = 2Ti5{^ + Vp) (29) 
-oo 

which arises from the 0^0"+ term in ^j™"* Eq.(0) when acting on the motionally cooled ground 
state Eq.Q. For blue sideband detuning A = — |A| < we obtain a non-zero contribution 
rom the delta function. In the usual ion trap excitation schemes for quantum computing 
13^ . the resonant (rotating wave approximation) portion of the Hamiltonian Eq.(jS)) gives 
rise to the anti-Jaynes-Cummings type interaction Eq.(jS)), which permits transitions of the 
form \n)\g) — |n + l)|e). With our motionally cooled ground state Eq.ijH} with ri = such 
transitions are possible from \Q)\g) Thus ap(7+ represents a resonant contribution 

to the Hamiltonian. 

For red sideband detuning, A > we do not obtain a contribution from the delta function. 
Here the resonant portion of the Hamiltonian Eq.(0) gives rise to the usual Jaynes-Cummings 
type interaction Eq.((7j), which permits transitions of the form \n)\g) ^ \n — l)|e). With an 
initial motional ground state Eq.Q with n = such transitions are not possible. However, 
in writing down Eq.(|7j) we have dropped the non-resonant anti-normally ordered terms. It 
is the aj,cr+ non- RWA term in ^j"^ Eq.© that can produce transitions from the motionally 
cooled ground state Eq.(jni) and is responsible for the Unruh-like behavior of the probability 
Pm for excitation out of this ground state. 

We can also give an Unruh analogy interpretation of the zero contribution to the delta 
function 5(A + z/p) in Eq. lj^ . The case of red sideband detuning A > 0, a constant trap 
frequency is analogous to an inertial observer moving with constant velocity in Minkowski 
space, giving rise to a constant Lorentz transformation as discussed above. For a constant ve- 
locity observer, an inertial detector would have a response function proportional to 5{E + Eq) 
where E (analogously A) is the energy of the detected particle and Eq > Q (analogously Vp) 
represents the atomic transition frequency producing the particle to be detected [3|. Since 
the detector only responds to positive energies > 0, the contribution from the delta func- 
tion 5{E + Eq) is zero, which simply states that an inertial detector moving through the 
Minkowski vacuum would detect no particle production. 
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III. FINITE CHIRP, GENERAL EXPRESSION 



For a finite chirp between times to and T we can develop a general expression for 
Pm(T,to) in terms of the incomplete gamma functions 7(/i, x) = Jq dt t^~^e~* and 
r(/i, x) = fl^dt t'^^^e^* such that 7(/i, x) + r(/i, x) = r(/i). Let us write Jy^ dy = 
f^dy — Jq" dy — f^dy for general finite duration limits in Eq.()22|) with the definitions 
yo = e'^^° and yT = e'^'^. By scaling the integration variable to y = y^x in the second integral 
on the right hand side and y = yrx in the third integral on the right hand side one can 
easily show that 

1 Tiia] ( *iV'^ 

/p(T, to) = - ^'""^ f/' (l - I'ita, -ihy^) - V\ia, -ibyr)) , (30) 



with the same definitions of a and b used in Eq.(j22)), and where we have defined the nor- 
malized gamma functions 7'(/i, x) = 7(/i, x)/r(/i) and r'(/i, x) = r(yU, x)/r(/i) such that 
7'(/i, x) + r'(yU, x) = 1. Thus we obtain 



p=l ^ ' ^ ' 



2|^(P)|2 



(31) 

The previous expression for the total excitation probability in Eq.(|25|) corresponds to 
Pm{oo, — oo) using Eq. (j3ip above which formally corresponds to sweeping the trap frequency 
from an initial zero value to an infinite final value. Considering a more realistic situation 
applicable to experiments, let us consider to = and a finite trap expansion time T, which 
corresponds to sweeping the trap frequency from u —>■ ue'^'^. Considering Eq. ()31|) as a 
function of the detuning A with parameters Up/n and nT we can recover Eq.()25|) under the 
following conditions 

^«1, ^e^^»l Vp=(l,2,...,iV) (32) 

which makes the incomplete gamma functions small compared to unity. As an example, 
taking Up/n = 0.01 and {vp/ k) exp(KT) = 100 requires that kT > 9.2 

Note that Eq. (j^T|l approaches Pm{oo, — oo) of Eq. lj^ as to ^ ~oo and T — > oo in which 
the incomplete gamma functions approach zero. The main point for experimental purposes 
is that we only need to take finite limits of the order Up/n ~ 10^^ and {vp/ k) exp(KT) ~ 10^ 
to approximate the full Unruh case, Pm(oo, — oo) as shown in Fig.Q. 
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RT,(t ,-oo) 




X=27tA/K 



FIG. 3: Finite duration chirp: Plot of Eq. (|31|) for p = 1 and ^Iqij/k = 1, Up/n = 1 with x = 2ttA/k 
as the abscissa and Pm{t, — oo) as the ordinate for various values of yt = e'^'^ = vi{T)/u; red curve 
t/T = 1) green curve yx = 10, blue curve yx = 100. The smooth black curve is the pure Unruh case 
Pm{oo, — oo) given by Ea. ((^ . 

IV. DISCUSSION AND CONCLUSION. 



As we have shown, the experimental signature of the exponential modulation of the trap 
frequency is the Planck-like form for the excitation probability for the two level electronic 
system in each ion. In such experiments it is the ratio of the excitation probability on the 
red (A = v) and blue (|A| = z/) sidebands that is determined: 

Re = ^ (33) 

m 

as this number is independent of the Rabi frequency, the Lamb-Dicke parameter, and the 
time of interaction between the vibrational and electronic degrees of freedom IJ, |l5| . Let us 
consider the case of a single ion with trap frequency u. Using Eqs. ()25|28|l we see that 



R 



(34) 



In a typical experiment one can detect R values as low as 0.05 with about 20% error. This 
implies that v / k, 0.5. At secular frequencies of v/2t[ = 0.1 — 1 MHz, we need a modulation 
frequency of the order of a few hundred Khz to MHz; a not particularly difficult requirement 
for fast electronics. 
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The key issue however is the absolute size of the excitation probabihties at the red and 
blue sideband. This is determined by the prefactor {Qoi])'^ / [ku) . Defining 

z = 2711'/ K, (35) 

the equations for the excitation probability at the red and blue sidebands are 

(37) 



As we expect z to be of the order of unity, we require that the secular frequency is within 
one order of magnitude of the effective Rabi frequency. This corresponds to a rather weakly 
bound ion, but should be achievable if stimulated Raman transitions are used to couple the 
two level system. For example the relevant transition in 9Be^ can have an effective Rabi 
frequency of the order of Qo/2n = 500kHz If we use the centre of mass mode with 
secular frequency of z//27r = 200kHz, and a Lamb-Dicke parameter of rj = 0.2, the prefactor 
is 0.25. At a more conservative estimate of rj = 0.05, well within the Lamb-Dicke regime, 
the prefactor drops to a value of 0.015. These numbers are encouraging enough to suggest 
the plausibility of observing an analogous Unruh-like effect in today's linear ion traps. 
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